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TECHNICAL NOTE 2597 

INVESTIGATION OF LAMINAR BOUNDARY LAYER IN COMPRESSIBLE 
FLUIDS USING THE CROCCO METHOD 
By E. R. Van Driest 


SUMMARY 


In the present investigation of the flow of air in a thin laminar 
boundary layer on a flat plate, the Crocco method has been used to solve 
the simultaneous differential equations of momentum and energy involved 
in such flow. The Crocco method was used because it gave accurate 
results for arbitrary Prandtl number near unity. The Prandtl number 
was taken at 0.75, the specific heat was held constant, and the Sutherland 
law of viscosity- temperature variation was assumed to represent the Q 
viscosity data starting with an initial ambient temperature of -67-6 F. 
The main results presented here are the skin-friction and heat-transfer 
coefficients as functions of Reynolds number, Mach number, and wall-to- 
free-stream temperature ratio. Variations of shear, velocity, tempera- 
ture, and Mach number across the boundary layer are included. The 
Crocco method is discussed in detail. 


INTRODUCTION 


In the present investigation of the flow of air in a thin laminar 
boundary layer on a flat plate, the Crocco method (reference l) has been 
used to solve the simultaneous differential equations of momentum and 
energy involved in such flow. The Crocco method was used because it 
gave accurate results for arbitrary Prandtl number near unity. The first 
assumptions for the present discussion are as follows: 

( 1) The boundary layer is thin 

(2) The pressure gradient is zero 

( 3 ) The plate surface is smooth 

( 4 ) The motion is two-dimensional and steady 

( 5 ) The Prandtl number is constant 
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SYMBOLS 


Variables, parameters, and functions 

x,y geometrical coordinates parallel and perpendicular to flat 

plate, respectively 

u,v velocity components of boundary- layer flow in x- and 

y-directions, respectively 

p mass density 

[i coefficient of viscosity 

k coefficient of thermal conductivity 

T absolute temperature 

i enthalpy 

Cp specific heat at constant pressure 

7 ratio of specific heat at constant pressure to specific 

heat at constant volume 

t shear stress 

M Mach number 

Pr Prandtl number ^Cpp/k^ 


Poo^ 

a) exponent of viscosity- temperature power law 

S Sutherland constant 

5 thickness of boundary layer 


shear function 
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i,j je 1 ,© 11 

5 e* 

Roo 

c f 

c f 

q 

h 

C H 

r 

** = i / i oo 
% = 

p » * p / p « 

= ^/^oo 

M* = M/M to 
T* = T/T^ 


ratio of Sutherland constant to free-stream temperature 

(no 0 k/tJ 3 k) 


functions connected with enthalpy distribution 

shear function near outer edge of boundary layer 

correction to near outer edge of boundary layer ( in 

computations, at u* = O. 98 ) 

Reynolds number (p^x/^ 

local skin-friction coefficient (2 t w fcu^ 2 J 


mean skin-friction coefficient 



T v /Poo U oo X ' 


rate of heat transfer 

local heat-transfer coefficient, difference between u* and 
unity as u^ approaches unity 

Stanton number 

recovery factor 


distance parameter 



f(u*) = 2u*p*p* 
g 2 (u) = t \[2x 
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Subscripts 


oo free-stream conditions 

w wall conditions 


TRANSFORMATION OF BOUNDARY-LAYER EQUATIONS 


The equations governing the steady two-dimensional motion of a 
compressible fluid in a thin boundary layer on a flat plate with zero 
pressure gradient are: 


S (pu) +|;(pv) = 0 


( continuity) 


( 1 ) 


Su 

pU ^ 


du 



(momentum) 


( 2 ) 


pu 




( energy) 


(3) 


Because of the thinness of the layer, Sp/Sy = 0, and hence for a per- 
fect gas the density varies inversely with the temperature. Upon 
assumption that the Prandtl number Cp|a/k is constant and use of 

0 ^ = di/dT, equation (3) becomes 



The independent variables are first transformed from x and y 
to x and u by means of 


u = u(x,y) 


Use of t = p Bu/dy and the elimination of pv then yield the fol- 
lowing two equations: 
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u 


d / PH \ 
dx\ t ) 



(5) 


(1 - Pr) 3s T ^ + (§ + Pr ) r2 • pt, »r 0 (6 > 

with boundary conditions 



T = 

00 

at 

0 

ii 

X! 

0 < u < 

(7a) 

St 

O 

II 

> 

•H 

II 

•H 




(7b) 


at 

X > 0 , 

u = 0 

T 

= 0, 

i = i 

00 

at 

X > 0 , 

u = u 

00 

(7c) 


It is now assumed that the enthalpy is a function of u only; that is, 
di/dx = 0. (This can be shown (reference 2) to be the case for a 
Prandtl number equal to unity, whether the wall is insulated or not.) 
Hence it follows from equation (6) that a type of solution is 
t = g^(x)gg(u) . Whence, since p and |i are functions of temperature, 

and therefore enthalpy, and therefore velocity u alone, there results 
from equation ( 5 ) and boundary condition (7a) that t = g2(u)/\/2x. 
Equations (5) and (6) now become 

g 2 g " 2 + U PP = 0 (8) 


(i" + Pr)g 2 + (l - Pr) i'g’ 2 = 0 (9) 

with boundary conditions 

g' 2 = 0, i = at u = 0 (10a) 

g 0 = 0, i = i at u = u (10b) 

d 00 00 


where the prime signifies differentiation with respect to u. When the 
variables i, u, p, and p are made dimensionless by putting 

0 f u * = u / u °°-> P* = P/Pco> and- P* = P/Po°> is seen that S* 


d* = i / i o 


must then be equal to 


p ll u 

*oo f oo 00 


- gg(a) . Equations (8) and (9) become 


finally 
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g*g*” + 2u*p^ # = 0 


( 11 ) 



( 12 ) 


with boundary conditions 


S*' = °> i* = ^00 at u * = 0 


(13a) 


g* = °, i* = 1 


at u. 


1 


(13b) 


The factor 2 is inserted in equation ( ll) so that when = 1> equa- 

tion (ll) will become the Blasius equation (reference 3) > the exact 
solution of which will be used as the first step in the Crocco method 
of solving equations (ll) and (12) in the case of gases. 

The boundary- layer equations (equations (l), (2), and (3)) have thus 
been transformed into two simultaneous ordinary differential equations 
which lend themselves more readily to solution. Having the solution of 
equations (ll) and (12), one can then compute the coefficients of skin 
friction and heat transfer and also the variation of certain properties 
across the boundary layer. 


SOLUTION OF EQUATIONS FOR ARBITRARY PRANDTL NUMBER NEAR UNITY 


In order to solve equations (ll) and (12), it is necessary to 
assume a law of variation of viscosity with temperature. Two laws 
have been frequently used to represent viscosity data, namely, the 
power law and the Sutherland law. The power law is 


where, in the case of air, o> ranges from O.76 for ordinary temperatures 
to 0.5 for higher temperatures. The Sutherland law is 


AND ARBITRARY VISCOSITY LAW 



(l4) 




T = T/T 

* ' 00 
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= T 1 / 2.-1 + 6 

* i + (e/n*) 


6 = S/T 


(15) 


where S is a constant depending upon the gas used. From data of 
Tribus and Boelter (reference 4) , the constant S is found to be 110° K 
for air. Both laws are plotted in figure 1 using initial conditions 
corresponding to the standard isothermal altitude range of 35,332 
to 104,987 feet. It is seen that the Sutherland law fits the data much 
better than the power law; therefore, since both o> and Q depend 
upon the free- stream temperature T^, it follows that the Sutherland 
law is most suitable for use in equation (ll). When it is assumed that 
the specific heat is constant and i^ is set equal to CpT^, then i* 

can replace T* in equation (l4) or (15), and p* = i*" 1 . If the varia- 
tion of the specific heat is to be taken into account, then it is only 
necessary that the relation between enthalpy and temperature (from gas 
tables) be available in order to convert the viscosity and density to 
functions of enthalpy. 

Following Crocco, equation (12) is written in the form 

^-(VK» 1 ' Pr ) + Pr f 2 - e* 1 " Pr - ° 


the solution of which is, in terms of and boundary conditions i*(0) 


and 


'(0), 


i*(%) = i*(0) + V(0)l(u*) - Pr^- J( % ) 


(16) 


where 


I = 


r* 


lo 



-.Pr-l 


du„ 


J = 


r* 

r s* _ 

>0 

Le»(o)_ 


-Pr-l 


chr. 


r*i 

r f 

1 0 

k<°>J 


-1-Pr 


du„ 
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This solution can be written in terms of the boundary condition in 
equation ( 13b) instead of i*'(0) by letting i* = 1 at u* = 1 in 

equation (l6) and eliminating i*'(0) between the resulting equation 
and equation ( l6) . Thus 


u„ 


1 = i (0) + i '(0)1(1) - Pr -A J(l) 

1 m 


( l6a) 


from which and equation ( l6) , 


i*(u*) = i*(0) - 


I(u») 


i*(0) - 1 + Pr - — 


u oo 2 

p(%) 

kT 

Ll(l) 


J( 1) - J(u*) 


or 


. , V 1 W / 1 W I.L V U oo _II/ \ 

i*(%) = -r- - (— - IJB (uj + — 0 (%) 

- L oo \ x oo / - L 00 


(17) 


where 


0 I (u*) 


- I ( u ») 

"KlT 


6 II (u x .) = Pr 


0 I (u*)J(1) 


J(%) 


The method of successive approximations suggests itself for calcu- 
lation of the shear function from equations (ll) and (17)* That is, a 
first approximation to g* would yield from equation (17) an i* which 

in turn would lead to a new through solution of equation (ll), upon 

satisfaction, of course, of the boundary conditions in equation (13)« 

The new g* could then be used to compute another i* which would give 
another g^, and so on. 


After extensive calculations Crocco observed, however, that the 
enthalpy distribution i*(u*) was practically independent of the law 
of variation of with enthalpy, that is, essentially independent 

of the law of variation of viscosity with temperature . This can be 
deduced upon inspection of the expressions for I and J in which g* 
and g*(0) appear as a ratio. It is apparent that although a change 
in the viscosity-temperature law would affect the absolute value of g* 
in equation ( ll) , nevertheless such a change could hardly affect the 
value of g*^(0), because g*(0) would be affected in much the same 
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way as g* because of the variation in viscosity law. Furthermore, for 
Prandtl numbers near unity, the effect of viscosity- law variation on 
enthalpy distribution would be even less. In his numerical proof of 


values of 0 from 0 to 3 when i # = in equations ( l4) and ( 15 ) ■ 

As a result of the above observation, one concludes that it is not 
necessary to carry through a lengthy iteration process between equa- 
tions (ll) and (17) • Rather, the enthalpy distribution can be immediately 
computed using the shear function corresponding to any viscosity law. 

The final shear function is then obtained as the solution of equa- 
tion (ll) by using an exact viscosity- temperature law. This simplifica- 
tion of solution of the simultaneous equations (ll) and (17) is the crux 
of the Crocco method of laminar -boundary-layer analysis. 

The shear function which may be readily used to compute the 
enthalpy distribution is that corresponding to the viscosity law such 
that P*h* = 1. (Since p* = T*"l, it is then necessary that p* = T*, 
whence o = 1 in equations (l4).) Therefore, equation (ll) reduces to 
the Blasius equation which has already been solved exactly. 

Table I gives the Blasius solution of equation (ll) when P*d* = 1» 
Table II gives the values of I and J based upon the Blasius solution; 
table III gives the values of 0^ and 0^. 

It may be well to interject here that the special case P*P* = 1 
represents a gas or state of a gas for which there is no effect of com- 
pressibility on skin friction for thin laminar boundary layers, although 
the temperature of the fluid rises because of dissipation of the kinetic 
energy, and properties of the fluid vary accordingly. 

Now that the enthalpy distribution has been calculated, attention 
will be directed toward the details of solution of equation (ll) • 

Writing equation (ll) in integral form, taking into consideration the 
boundary conditions ( 13 ), one obtains the integral equation 


-in which p p is a known function of u depending upon the viscosity 
* * 

law assumed. The method of successive approximations will again be 
utilized. However, as Crocco pointed out, the direct method of 
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substitution, integration, resubstitution, and reintegration does not 
converge upon a single solution but yields values of g* which oscillate 
about the exact value. For, if an initial value of the shear func- 
tion equal to Ag* , in which A is a constant and g* is the 

exact value, is substituted into the right-hand side of equation (l8), 
the new g^ obtained upon integration will be equal to g„ /A, resub- 


*2 


S/ As 


stitution of which will yield = Ag = g Therefore, it is seen 

*3 *e *1 

that the next substitution should be \/g* g* = g* • In the iterative 

V 1 2 e 

process of solving equation (ll) (equation (l8)), it then follows that 
successive values g^ for resubstitution into equation (l8) should be 
the geometric mean of the two previous substitutions. Although an 
exact viscosity-temperature relation is now used in equation (l8), the 
iterative process may be started with the Blasius solution (table I) 
which corresponds to a linear viscosity- temperature law. 


It is next observed that a singularity exists at u* = 1 owing to 
the boundary condition (l3b) that = 0 at that point. Therefore, 
integration is carried out only to an upper value of u^ = 1 - h, 
where h is an arbitrarily small value greater than zero. Equation (l8) 
becomes then 


g*(u*) 


g*(l - h) 


+ 



f(u*) 

>o e*( u *) U * 2 


(19) 


where f(u*) = 2u^p^n^. In order to avoid double integration, equa- 
tion ( 19) is written in the form 


g*(u*) = g*(l - h) + (l - 


U u 



f (u*) 


du * 


h 



f(u*) 

S*( u *) 


du v 


+ 



(l - u*) 


f(u # ) 

W iu * 


( 20 ) 


Since the first term g*(l - h) on the right-hand side of equa- 
tion (20) is unchanged by successive iteration, it is necessary to 
derive a method to adjust g*(l - b) in each iteration so that the 
boundary condition g^(l) = 0 is more nearly approached. Hence the 
problem is now reduced to that of determining an accurate value 
of g*(l - h) which is compatible with the boundary, condition g*(l) = 0. 
A method of estimating g*(l - b) is now given. 
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Following Crocco, one notes that for u* — > 1 equation (ll) 
approaches 


g * g * 


+ 2 = 0 


( 21 ) 


which can he integrated in closed form. The double bar is used to 
designate the solution of this approximate equation. The first integral 
of equation (2l) yields 


g^' = -2\/log e (B/g*) 


( 22 ) 


where B is an arbitrary constant and the minus sign is chosen because 
of the sense of g^' . Further integration gives 


1 





V 10 ® J B /eP 


(23) 


in which the boundary condition g*(l) =0 is satisfied. Upon substi- 
tution of log (B/gp = \r , this equation can be expressed in terms of 

2 r x -t 2 

the tabulated error function 0(x) = -j= / e dt, thus 

V Wo 


1 


u„ 


1 

2 


B 


\/rf Jl 





(24) 


The constant B can be eliminated between equations (22) and (24). 
Hence 


(l - u*)/g„ - 4^ e 


(I^72) 2r 


1 - 0(-|“V2) 


( 25 ) 


from which S*(l - h) can be calculated, given h = 1 - u* and the 
slope g^'(l " h) • Since the slopes of and are to match 

at u = 1 - h, one puts 
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g ' (l - h) = g * ( 1 - h) ( 26 ) 

* * 

which is computed from the first integral of equation (ll), namely, 


g*'(l - h) 


/ 1-k f(u ) 

J 0 du * 


(27) 


The ratio h/g^(l - h) is plotted against g*'(l - h) in fig- 

ST(1 ~ h) 


ure 2. Crocco noticed that the combination , 


- 1 


Ig^'U “ h)| 

can be represented by the straight line O.7828 + 0.0178 |g^'(l - h) | 
over the practical range of g^'(l ~ h) between -2.2 and -4 (fig. 2). 

Now the error between |^(l - h) and g (l - h) will depend upon 
the smallness of the h, that is, upon how well the approximate equa- 
tion (21) represents the exact equation (ll) . As Crocco points out, i*, 
and therefore varies rapidly for Pr < 1 when u* approaches 

unity (cf. tables II and III and equation (l7))j hence it may be 
necessary to apply a correction so that 


g (l - h) = g (1 - h) + 5g 


(28) 


Rewriting equation (l8) in the form 


K.fu.) = (1 - f * (1 - iylj du * (29) 

* 


and combining with the integral of equation (ll), namely, 


g '(u ) = 
* * 


r% au . « • 

0 g*( u *) * * 


(1 - h) 


f(u*) 

/!.„ du ’ 


there is obtained 
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r u* f(uj 

Sx^d*) = -(l - U^g^'Cl - h) + (l - u*) / / \ du# + 

J 1-h s * y *' 



f(u # ) 

(1 - -> 


du# 


Similarly from equation (2l) one obtains 


(3D 


g#(%) = -(1 - u*)g*'(l - h) + (1 


U„ 


uJ / 2 


1-h ^*( u *) 


du * + 




du * 


Subtraction of this equation from the preceding one gives 


(32) 


&g*(u*) = (l - u # ) 




(1 - 


u 


. /- - 
*\e» 



(33) 


which can be solved in principle by iteration since f is a known func- 
tion of u*. That is, starting with an assumed g* , g*^ is avail- 

able from equation (25), whence 5g#.^(u*) is computed, whence 
g* 2 (%) = d*) + Sg* 2 (u*), and so forth. 

However, the above iterative process can be rearranged as follows: 
Write equation ( 33) thus 


5 s*(%) 




f 2 | f \ + 

~p r ~ IF” 7F~ I 
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■nU. 


Sg*(u*) = (l - u*) 


* f - 2 
'l-h ^ 


du* + / (l - u*) - 2 du* + 

- s; 


■Al . 


(1 - u*) / * hr -iV d % + / (1 - %)hr -=i f 


l-h I s * g. 


u 


g* 


The term — 
8* 


1 1 


- = » can now he expanded in the form 
i i &g* & s* 2 

«* ""f= 2 r 3 " ' ' * 


&&* 


5 S* 


whereupon, when -'jt is neglected since * =r is small, equation 

«; 3 

becomes 


&g*r(u*) = (1 - u*) / * 2 du* + / (l - u*) f _I 2 du* 


l-h R 


u 


P u * 5g pi 5g 

(1 - u*) / -^zr| f du* - / (1 - U*) f du* 


'l-h 


u* 


£* 


If the first approximation to Sg* is obtained by letting Sg* = 0 
the right-hand side of equation (37)> then 


r * f _ 2 ru f - 2 

~ =s= — du* + / (l - u*) -= — du* 

h S* ^u* S* 


u * 8g* i 


6 g* n 


6g*p(u*) = Sg* (U*) - (l - U*) / - m r T f du * ~ (l “ M*) ~ 

2 1 Jl-hF d ^ u„ g.. 


(35) 


(36) 


(35) 


(37) 

on 


f du* 
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&S* 3 (u*) = Sg^u*) - (1 - u*) 


'% 5g *c 


ol 


1 - h ^ 2 


f du* - 


(1 - u*) 


% 


5g * £ 

=2 

S * 


f du* 


and so forth. Stopping with the third approximation, the value 
of 5g^ at u* = 1 - h is then 







5 S*p 

(l - u*) — f du # 
g * 2 


(38) 


Crocco now gives an approximate expression for equation ( 38 ) in 
the form 




&g*(l - h) = 2~<a(i # - l) 


S„ 


4 et( 


> - "'If • MB' 


(if ♦ • • • 

+ h(i* - l) 2 

/h \ 2 

B o - B lM= + • • • 

Wl 


L W / J 


:(i* - l)- 


c 0 - c 7 i \ 2 + . . . 


- h 


B o - Di w 


where 


(39) 


1 1 
a = — - 


h = - 


2 0+1 
1 a_ 

B " 0 + 1 


1 h 

C IB ‘ 0 + 1 

are the coefficients in a Taylor series expansion of about the 

point i = 1 using the Sutherland law (equation (15)) for 
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viscosity- temperature variation. The other coefficients in equation (39) 
are 


A„ = 


o 3 - (2 - Pr)T! 


A-, = 


A 0 = 


2 

1 

2 - (2 - Pr)t) 

_1 + 2(1 - T] ) 


4 

S’ 

1 

2 - (2 - Pr)V 

jl + 2(1 - T| ) 

^0 

1 + 4(1 - Tl) 


5 - (4 - Pr)ii 


7 - (6 - Pr)Ti 


An - 


3 ~ 2 - (2 - Pr)Ti 


1 _ A o 

4 - (3 - Pr) q 6 - (5 - 2Pr)^ 


B o 4 - (3 - 2Pr)t) 


B n = 


r i 


B, 


1 3 - (3 - 2Pr) r| [l + 2(1 - q) ' 6 - (5 - 2Pr)n 


= 


0 " 5 - (4 - 3Pr)T) 


Ci = 


2 

1 

4 - (4 - 3Pr)ri 

[l + 2(1 - ti) 


B, 


= 


0 3 - n 

2 


D-, = 


1 2 - T\ 


1 + 2(1 - T ) ) ' 2 + 3(1 - 4) 


Vi — - ___ 

where ti =-— 17 ' • In equation (39)> the values of g*, g^', and i # 
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are to be determined at u* = 1 - h. It will be found that the terms 
included in equation (39) will give sufficient accuracy to Sg^(l - h) . 
In fact, in the calculations carried out in the present report it was 

found that 6g* was negligible, being of the order of 10“^. 


SKIN-FRICTION AND HEAT-TRANSFER COEFFICIENTS 


62 ( u ) 

Since r = and g* 

\j£x 

shear stress t and the shear 


p Ll u - 
_ 00 00 00 

function 


ggCu), 
g* are 


it follows that the 
related by 


r 



where R^ = — - °° , the Reynolds number in terms of the length x from 

^00 

the leading edge to the point in question on the plate. The surface 
stress is then 


T l pu 2«*<£i 

w 2 ““ 

In terms of the local skin-friction coefficient Cf defined by 


(4o) 


T 


W 


c 


f 


2 

00 


(4i) 


there results 


c f ^=e*(o) (42) 

Furthermore, in terms of the mean skin-friction coefficient Cf defined 
by 



C 


1 2 
f 2 p oo u °o x 


(43) 




( 44 ) 


there follows 
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Comparison of equations (42) and (44) shows that regardless of Mach num- 
ber, Reynolds number, or wall-to-free-stream temperature ratio 

C f = 2c f (45) 

The heat transferred from the wall to the fluid per unit width of 
plate is given by 


M, 

V(0) T w 



Ww* 

k w 

loo 

C p dw U oo 

w 


iJ 


Pr u * 
00 


1 (0) T w 


From equation (l6a) one has 


V<°> = IHT 


1 - i # (0) + Pr J(l)^2- 


so that equation (46) becomes 


<3w = 


Pr 


2 T [ 

KD W L 


u 


-a r w 1 + Pr J(l) -^2- - i*(0) 


for variable specific heat, or 


q — - t -<T 

w Pr 1(1) w 1 “ 




1 + 2Pr j(1)Z-LJz M 2 


- T 


w 


(46) 


(47) 


(48) 


for constant specific heat. 
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Defining the local heat-transfer coefficient by 

1 


h = 


c Poo T v 


Pr 1(1) u* 


(^ 9 ) 


and the Stanton number by C 0 = h/(^ c p OT p oo u oo) > there follows, utilizing 
equation ( ll) , that 


C H “ 


c f 


Pr 1(1) 2 


(50) 


It is now interesting to observe from table II, that, closely, 


Pr 1(1) * Pr' 


2/3 


and 


2Pr J(l) * Pr 


1/2 


so that practically 


C H = 


Pr 


273 2 


and 

q w = -ChP^u 
when the specific heat is variable, or 



1 + Pr 1 / 2 7 ^ 1 M 2 ) - T 

2 00 j w 

when the specific heat is constant. (Specifically, for 
Pr 1(1) = 0.8302 and Pr 2 / 3 = 0.8258; 2Pr j(l) = 0.865^ 
and Pr"*"^ 2 = 0 . 8660 .) 


"*w 


= -0 w c p u 

n p 00 00 


T *( 


( 51 ) 


( 52 ) 


(53) 


Pr = 0.75, 
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Suppose now there is no heat transfer into or out of the boundary 
layer ( insulated wall) . Then q w = 0 and 

/ u 2 

i w = i^ + Pr 1 ' 2 (variable c p ) (5*0 

or 

T w = T^l + Pr 1//2 7 ~ 1 M m 2 j (constant c p ) (55) 

which is the wall temperature due to viscous action alone. However, 
the "recovery factor" is defined hy 


r 



(56) 


Hence, it follows upon comparison of equations (5*0 and ( 56 ) that, 
closely, 


Recovery factor = r = Pr 


1/2 


(57) 


for laminar boundary layers. 

In -the determination of the friction and heat transfer at the 
plate, it was not necessary to know the internal geometric distribution 
of properties of the boundary layer. However, if the geometric distri- 
bution of properties, such as velocity and temperature in a stability 
analysis, are desired, then a suitable distance parameter is 


since 




(58) 


RESULTS OF CALCULATIONS 


Calculations of the skin-friction and heat-transfer coefficients 
for laminar boundary layers on flat plates have been carried out for 
Mach numbers up to 20 and wall- to- free -stream temperature ratios of l/4 
to 6. Although the results for the higher Mach numbers are academic 
because of the assumptions made concerning the variation of properties 
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of the gas at high temperatures, these calculations were performed in 
order to indicate trends and compare with more accurate results. The 
Prandtl number was taken at 0.75, the specific heat was held constant, 
and the Sutherland law for viscosity-temperature variation was assumed 
to represent the viscosity data starting with an initial ambient tempera- 
ture of -67 . 6 ° F. The ratio of specific heats was taken at 1.400. 

The calculation procedure was as follows: The enthalpy distribu- 

tion was calculated using the Blasius distribution of the dimensionless 

shear function as given in table I; that is, values of 6 ^ and ©H from 
table III were used in equation (17) • The final shear distribution was 
computed by integration of equation (ll) by the method of successive 
approximations from u* = 0 to u* = O .98 starting with the Blasius 
solution in table I. In the iterative process, the value of the shear 
function at u* = O .98 was obtained from the solution of equation (21) 
(fig. 2). The required accuracy on g*(u*) was within 0.0025, hence 
within 0.5 percent for friction and heat transfer at the plate over the 
variation of wall- to-free- stream temperature ratios chosen herein. It 
was unnecessary to apply the correction 5g* of equation ( 28 ) because 

6 g^ was always found to be very small (of the order of 10"^). A work- 
sheet for carrying out an iteration for g* is given in the appendix. 

In figure 3 is plotted the mean skin-friction coefficient as func- 
tion of free-stream Reynolds number, free-stream Mach number, and wall- 
to -free- stream temperature ratio. Figure 4 gives the local heat- 
transfer coefficient as function of the same dimensionless variables. 

In figure 5 are compared results obtained using the Sutherland and power 
viscosity laws. Figure 6 shows the variation of the thickness of the 
laminar layer as a function of the above independent variables; the 
thickness was arbitrarily taken at the velocity ratio u* = 0.995 
obtained by estimating the g* distribution to that point. 

Figures 7 to 12 show the dimensionless shear-function distribution 
for various Mach numbers and temperature ratios including the insulated- 
plate case. Figures 13 to 15 give the variations of the velocity 
ratio u*, the temperature ratio t/Tqo, and the Mach number ratio m/m^ 

as functions of free-stream Mach number M^ and distance parameter 

\/Roo in the case of the insulated plate; figures 1 6 to 30 show the 

same dependent variables as functions of the same independent vari- 
ables in the case of wall- to-free- stream temperature ratios l/4, 1 , 

2, 4, and 6 . Figures 31 to 47 show the same dependent variables as 
above but as lines of constant temperature ratio instead of as lines 
of constant free-stream Mach number, the distance parameter now being 
the relative distance y/b. For purpose of comparison with figures 13 
to 15 the internal characteristics of the boundary layer are plotted 
in figures 48 to 5 0 tor Pr = 1 . 
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Figure 51 shows the variation across the boundary layer of the 


ratio of the total enthalpy per unit mass 


( 



at any point 


to that of the free stream for an insulated plate. For Pr = 1 the 
ratio would he constant at unity. The fact that Pr is less than 
unity in these calculations brings about the transfer of the energy 
from the inner to the outer layers of the boundary layer; this is due 
to the readjustment of the relative effects of viscosity, heat capacity, 
and heat conductivity in favor of the heat conductivity. 


Aerophysics Laboratory 

North American Aviation, Inc. 

Downey, Calif., January 9 > 1951 
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APPENDIX A 

WORK SHEET FOR CARRYING OUT AN ITERATION FOR g # 

The following headings and equations are used in carrying out an iteration for g*: 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

% 

o 1 

(table III) 

e 11 

(table III) 

i*(°) = £ 

u« 2 

i 

00 

i»(u,) 

(equation ( 17 )) 

U*(u*) 

(equation ( 15 ) 
and gas tables) 

P*0* 

(columns 6 and 7 
and gas tables) 

2 P*l 1 » u * - f K) 

(assumed) 


) 2 

^(u*) + T 2 - (column 6 ) 

1 co 


„ ( U ) = T i/ 2 - hm 


i + 0^05 


(column 7 ) 



U C o 

With constant specific heat, -2- = (7 - 1)M„Y and i* = T* for i*, = CpT^. 

2-CO 


11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

f(u») 

/"* T& du * 
Jq 

D 

( see below) 

E 

(see below) 

D - E 

H* H 

3 * 

1 

(H 

F 

(see below) 

g»(l - h) 

e» 2 K) 

(equation ( 20 )) 


8*3 = ^ 18*2 



. P u * f( % ) , 

= (1 -^X 4u * 


p!- h f(u#) 

Jo 


VII 


du * 


f K> 


(1 • u * ) g^T%T ^ 

g* g (u*) = g*(l - h) + D- E+ F 

g*(l - h) - ^(1 - h) + 6s» (Bg* is negligible in present calculations) 

1 r P* ( V h> - l] - 0.7828 + 0.0178 If;' (l - h)| 

|g7’(i - h)|L h J l*t 


— r 1 - 11 f 

*•<!-*> -v< l -h> r 


du 


2k 
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TABLE I 

Shear Function g u When p u = 1 

/V "/V 


u * 

e* 

u * 

e* 

u * 

g * 

0 

0.66411 

0.45 

0.61772 

0.875 

0.27994 

.05 

.66405 

.50 

.60013 

.900 

.23881 

.10 

.66361 

.55 

.57836 

•925 

.19293 

.15 

.66242 

.60 

.55183 

• 95 

.14097 

.20 

. 66009 

.65 

.51985 

• 96 

.11797 

.25 

.65625 

• 70 

.48161 

• 97 

.0933^ 

•30 

.65050 

.75 

.43607 

• 98 

.06659 

• 35 

.64245 

.80 

.38189 

• 99 

.03681 

.to 

.63167 

.85 

.31715 

1.00 

0 


ro 

ON 


TABLE II 


VALUES OF I AND J BASED UPON BLASIUS SOLUTION FOB g # 


U* 

Pr = 0.5 

Pr = 0.725 

Pr = 0.75 

Pr = 1.00 

Pr = 1.25 

Pr = 1.50 

Pr = 

2.00 

I 

J 

I 

J 

I 

J 

I 

J 

I 

J 

I 

J 

I 

J 

0 

.1 

.2 

• 3 
.4 

• 5 
.6 
.7 
.8 
.9 

• 92 
. 9 ^ 

• 96 

• 98 

• 99 
1.00 

0 

.1000 

.2002 

.3008 

.4025 

.5062 

.6135 

.7267 

.8505 

.9967 

1.0390 

I.0696 

1.1131 

1.1671 

0 

.0050 

.0200 

.0452 

.0806 

.1269 

.1849 

.2567 

.3454 

.4612 

.4899 

0 

.1000 
.2001 
• 3004 
.4014 
.5034 
.6073 
.7144 
.8268 

.9500 

.9770 

1.0056 

1.0363 

1.0708 

0 

.0050 

.0200 

.0451 

.0803 

.1260 

.1827 

.2512 

.3337 

.4343 

.4579 

0 

.1000 

.2001 

.3004 

.4012 

• 5031 
.6066 
.7130 
.8243 
.9451 

• 9715 

• 9991 
1.0286 
l. 06 i 4 

0 

.0050 
.0200 
.0451 
.0803 
.1260 
.1824 
.2507 
• 3324 
.4315 
.4542 

0 

.1000 
.2000 
.3000 
.4000 
.5000 
.6000 
.7000 
.8000 
.9000 
.9200 
.9400 
.9600 
.9800 
• 9900 
1.0000 

0 

.0050 

.0200 

.0450 

.0800 

.1250 

.1800 

.2450 

.3200 

.4050 

.4232 

0 

.1000 

.1999 

.2996 

.3988 

.4969 

.5935 

.6875 

.7776 

.8605 

.8757 

.8901 

.9037 

•9159 

0 

.0050 

.0200 

.0449 

• 0797 
.1240 
.1776 
.2396 
.3086 
.3817 
.3965 

0 

.1000 

.1998 

.2992 

.3975 

.4939 

.5872 

.6756 

• 7564 
.8251 
.8367 

.8472 

.8564 

.8638 

0 

.0050 
.0200 
.0448 
.0794 
.1232 
• 1753 
.2342 
.2980 
.3613 
.3731 

0 

.1000 

.1997 

.2985 

• 3951 
.4881 

• 5750 

.6531 

.7186 

.7660 

.7727 

.7780 

.7824 

.7852 

0 

.0050 

.0200 

.0447 

.0788 

.1213 

.1708 

.2246 

.2787 
.3276 
• 3353 

• 5591 

.5103 

.5050 

.4608 

.4245 

• 3941 

.3475 

.6369 

.7039 

.5602 

.5871 

.5523 

.5768 

.4901 

.5000 

.4429 

.4473 

.4060 

. 4 o 8 o 

• 3527 
• 3532 , 

1.2806 

1.1210 

I.1070 

.9252 

.8682 

.7863 
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TABLE III 


VALUES OF 6 1 


AND 6 11 BASED UPON BLASIUS SOLUTION FOR 


S* 



Pr = 

0.5 

Pr = 0.725 

Pr = 

0.75 

Pr = 

1.00 

Pr = 

1.25 

Pr = 

1.50 

Pr = 

2.00 

\1 

* 

e 1 

e 11 

e 1 

e 11 

e 1 

e 11 

e 1 

e 11 

e 1 

e 11 

e 1 

e 11 

e 1 

e 11 

0 

.1 

.2 

• 3 

.4 

• 5 
.6 

.7 

.8 

.9 

• 92 

• 94 

• 96 

• 98 
.99 

1.00 

0 

.0781 

.1562 

.23W 

• 31^3 
.3953 
.1*790 
.567^ 
.661*1 
.7783 

.8051* 

.8352 

.8692 

.9113 

• 9397 
l 

0 

.0250 

.01*50 

.0601 

.0703 

.0757 

.0761 

.0712 

.0610 

.01*31 

.0385 

0 

.0892 

.1785 

.2680 

.3580 

.1*1*91 

.51*18 

.6373 

.7378 

.81*71* 

.8715 

.8970 

.921*1* 

•9552 

•9733 

1 

0 

.031*1* 

.0615 

.0813 

.091*2 

.0998 

.0982 

.0891 

.0720 

.01*58 

.0390 

0 

.0903 

.1807 

.2713 

.3625 

.1*51*5 

.51*80 

.641*1 

.7446 

.8538 

.8776 

.9026 

.9292 

.9589 

.9760 

1 

0 

.0353 

.0633 

.0835 

.0966 

.1020 

.1002 

.0906 

.0728 

.0455 

.0390 

0 

.1000 

.2000 

.3000 

.4000 

.5000 

.6000 

.7000 

.8000 

.9000 

.9200 

.9400 

.9600 

.9800 

.9900 

1 

0 

.0450 

.0800 

.1050 

.1200 

.1250 

.1200 

.1050 

.0800 

.0450 

.0362 

0 

.1081 

.2161 

.3238 

.4310 

•5371 

.6415 

.7431 

.8405 

.9301 

.9465 

.9621 

.9768 

.9900 
• 9957 

1 

0. 

.0542 

.0958 

.1249 

.l4i4 

.1452 

.1367 

.1160 

.0842 

.0424 

.0336 

0 

.1152 

.2302 

.3446 

.4579 

.5689 

.6763 

.7781 
.8712 
.9503 
.9636 
• 9757 
.9863 
.9949 
.9981 

1 

0 

.0630 

.1109 

.1436 

.1611 

.1736 

.1509 

.1250 

.0862 

.0391 

.0301 

0 

.1271 

.2539 

.3796 

.5026 

.6206 

.7311 
.8305 
.9139 
.9740 
.9824 
.9894 
.991+8 
.9984 
• 999 4 

1 

0 

.0798 

.1394 

.1788 

.1974 

.1959 

.1749 

.1375 

.0882 

.0328 

.0234 

.0263 

.0235 

.0232 

.0192 

.0155 

.0124 

.0078 

.0123 

0 

.0081* 

0 

.0080 

0 

.0050 

0 

.0030 

0 

.0018 

0 

.0006 

0 


ro 
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Figure 1.- Viscosity-temperature data. 
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Figure 2.- Shear function g*(l - h) as a function of g^'(l - h). 
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Figure 3- - Mean skin-friction coefficient for laminar boundary layer of 
a compressible fluid flowing along a flat plate. Prandtl number, 0.75} 

e = 0.505. 


NACA TN 2597 






ii'iiliHHi! 


s aa a as 

Hi— MiBii 


Figure 4.- Local heat-transfer coefficient for laminar boundary layer of 
a compressible fluid flowing along a flat plate. Prandtl number, 0.75: 

e = 0.505. 
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Figure 5,- Comparison of mean skin- friction coefficients using Sutherland 

and power viscosity laws. 
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Figure 7 .- Shear distribution in laminar boundary layer of a compressible 
fluid flowing along an insulated flat plate at various free -stream 
Mach numbers. Prandtl number, 0.75> 0 = 0.505* 
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Figure 8.- Shear distribution in laminar boundary layer for wall-to-free- 
stream temperature ratio of l/4 and various free-stream Mach numbers. 
Prandtl number, 0.75j 0 = 0.505- 
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Figure 9*- Shear distribution in laminar boundary layer for vall-to-free- 
stream temperature ratio of 1 and various free-stream Mach numbers. 
Prandtl number, 0.75) 0 = 0-505- 
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Figure 10.- Shear distribution in laminar boundary layer for wall-to-free- 
stream temperature ratio of 2 and various free-stream Mach numbers. 
Prandtl number, 0.75j 9 = 0-505- 
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Figure 11.- Shear distribution in laminar boundary layer for wall-to-free- 
stream temperature ratio of 4 and various free-stream Mach numbers. 
Prandtl number, 0-75j 0 = 0.505* 
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Figure 12.- Shear distribution in laminar boundary layer for wall-to-free- 
stream temperature ratio of 6 and various free -stream Mach numbers. 
Prandtl number, 0-75j 9 = 0.505* 
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Figure 13.- Velocity distribution across laminar boundary layer on an 
insulated flat plate for various Mach numbers. Prandtl number, 0.75) 


e = 0.505. 
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Figure Temperature distribution across laminar boundary layer on an 

insulated flat plate for various Mach numbers. Prandtl number, 0.75$ 

e = 0.505. m 
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Figure 15- - Mach number distribution across laminar boundary layer on an 
insulated flat plate for various Mach numbers . Prandtl number, 0-75j 
6 = 0-505- 
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Figure 1 6 .- Velocity distribution across laminar boundary layer for wall- 
to-free-stream temperature ratio of l/k and various Mach numbers. 
Prandtl number, 0.75) 0 = 0.505* 
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Figure 1T-- Temperature distribution across laminar boundary layer for wall- 
to -free -stream temperature ratio of l/k and various Mach numbers. Prandtl 
number, 0.75) 0 = 0.505* 




NACA TN 2597 




Figure 18.- Mach number distribution across laminar boundary layer for 
wall-to-free-stream temperature ratio of l/k and various Mach numbers. 
Prandtl number, 0.75) 0 = 0.505* 
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Figure 19.- Velocity distribution across laminar boundary layer for wall- 
to-free- stream temperature ratio of 1 and various Mach numbers. Prandtl 
number, 0.75} 0 = 0.505. 
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Figure 20.- Temperature distribution across laminar boundary layer for 
wall-to-free-stream temperature ratio of 1 and various Mach numbers. 
Prandtl number, 0.75} 0 = 0.505* 
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Figure 21.- Mach number distribution across laminar boundary layer for 
wall-to-free -stream temperature ratio of 1 and various Mach numbers. 
Prandtl number, 0.75} © = 0.505* 
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Figure 22.- Velocity distribution across laminar boundary layer for wall- 
to-free-stream temperature ratio of 2 and various Mach numbers. Prandtl 
number, 0-75) 0 = 0.505. 
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Figure 23 .- Temperature distribution across laminar boundary layer for 
wall-to-free-stream temperature ratio of 2 and various Mach numbers. 
Prandtl number, 0.75) 0 = 0.505* 
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F Igure 24.- Mach number distribution across laminar boundary layer for 
wall-to-free -stream temperature ratio of 2 and various Mach numbers. 
Prandtl number, 0.75j 0 * 0.505* 






VJ1 


NACA TN 2597 










VJI 

ro 





Figure 25* - Velocity distribution across laminar boundary layer for wall- 
to-free-stream temperature ratio of k and various Mach numbers. Prandtl 
number, 0.75; 0 = 0.505* 
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Figure 26.- Temperature distribution across laminar boundary layer for 
wall-to -free -stream temperature ratio of b and various Mach numbers. 
Prandtl number, 0.75; 0 = 0.505. 
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Figure 27.- Mach number distribution across laminar boundary layer for 
wall-to-free -stream temperature ratio of 4 and various Mach numbers. 
Prandtl number, 0.75) 6 = 0.505* 
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Figure 28.- Velocity distribution across laminar boundary layer for vall- 
to-free-stream temperature ratio of 6 and various Mach numbers. Prandtl 
number, 0.75i 9 = 0.505. 
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Figure 29*- Temperature distribution across laminar boundary layer for 
wall-to-free- stream temperature ratio of 6 and various Mach numbers. 
Prandtl number, 0.75j 9 = 0.505* 
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Figure 30.- Mach number distribution across laminar boundary layer for 
wall-to-free- stream temperature ratio of 6 and various Mach numbers. 
Prandtl number, 0.75} 0 = 0.505* 
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Figure 31- - Velocity distribution a cross laminar boundary layer for = 0 
and various wall-to-free-stream temperature ratios. Prandtl number, 

0.75; 0 = 0.505. 
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Figure 32.- Temperature distribution across laminar boundary layer for 
= 0 and various wall-to-free-stream temperature ratios. Prandtl 
number, 0.75; 0 = 0.505. 
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Figure 33.- Velocity distribution across laminar boundary layer for 

= 4 and various wall-to-free-stream temperature ratios. Prandtl 

number, 0.75; 0 = 0.505* 
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Figure 34.- Temperature distribution across laminar boundary layer for 
= 4 and various wall-to-free-stream temperature ratios. Prandtl 
number, 0.75; 0 = 0.505. 
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Figure 35.- Mach number distribution across laminar boundary layer for 
= k and various wall-to-free-stream temperature ratios, Prandtl 

number, 0.75; 0 = 0.505. 
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Figure 36.- Velocity distribution across laminar boundary layer for 

= 8 and various wall-to -free -stream temperature ratios. Prandtl 
number, 0.75; 0 = 0.505. 
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Figure 37.- Temperature distribution across laminar boundary layer for 
M m = 8 and various wall-to-free-stream temperature ratios. Prandtl 
number, 0.75; 0 = 0.505. 
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Figure 38.- Mach number distribution across laminar boundary layer for 
= 8 and various wall-to-free-stream temperature ratios. Prandtl 
number, 0.75; 0 = 0.505. 


G\ 

VJI 


NACA TN 2597 







(T\ 

CT\ 



Figure 39*- Velocity distribution across laminar boundary layer for 

= 12 and various wall-to-free-stream temperature ratios. Prandtl 
number, 0.75; 0 = 0.505. 
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Figure 40.- Temperature distribution across 
= 12 and various wall-to-free-stream 
number, 0.75; 0 = 0.505. 
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Figure 4l.- Mach number distribution across laminar boundary layer for 

= 12 and various wall-to-free-stream temperature ratios. Prandtl 
number, 0.75; 0 = 0.505. 
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Figure lj-3.- Temperature distribution across laminar boundary layer for 

= 16 and various wall-to-free-stream temperature ratios. Prandtl 
number, 0.75; 6 - 0.505. 
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Figure 44.- Mach number distribution across laminar boundary layer for 
Mo, = 16 and various wall-to-free-stream temperature ratios. Prandtl 
number, 0.75; 6 = 0.505. 
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Figure 45.- Velocity distribution across laminar boundary layer for 

= 20 and various wall-to-free-stream temperature ratios. Prandtl 
number, 0.75; 0 = 0.505. 
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Figure k6.- Temperature distribution across laminar boundary layer for 

= 20 and various wall-to-free-strea.m temperature ratios-. Prandtl 
number, 0.75; 0 = 0.505. 
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Figure Vf.- Mach number distribution across laminar boundary layer for 
Mo, = 20 and various wall-to-free stream temperature ratios. Prandtl 

number, 0.75; # = 0.505. 
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Figure 48.- Velocity distribution across laminar boundary layer on an 
insulated flat plate for various Mach numbers. Prandtl number, 1.0} 
Q = 0.505. 
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Figure b9.- Temperature distribution across laminar boundary layer on an 
insulated flat plate for various Mach numbers . Prandtl number 1.0; 

e = 0.505. 
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Figure 50.- Mach number distribution across laminar boundary layer on an 
insulated flat plate for various Mach numbers. Prandtl number, 1.0$ 

e = 0.505. 
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Figure 51.- Total enthalpy variation across laminar boundary layer of a 
compressible fluid flowing along an insulated flat plate at various 
free-stream Mach numbers. Prandtl number, 0.75; © = 0.505. 


NACA TN 2597 






I 

I 


l 


